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MOLECULES  TO  THE  SECOND  VIRIAL  COEFFICIENT 

(REVISED) 

by 

Daniel  E.  Stogryn’and  Joseph  O.  Hirschfelder 
ABSTRACT 

"xi  * 

The  second  virial  coefficient  for  molecules  interacting  with  a 

spherically  symmetric  potential  is  divided  into  three  parts:  (1)  a  contribu¬ 
tion  B.,  related  to  the  equilibrium  constant  for  the  formation  of  bound 
b 

double  molecules;  (2)  a  contribution  B  re'ated  to  the  eouilibrium  con- 

m 

stant  for  the  formation  of  metastably  bound  double  molecules;  and  (3) 
a  contribution  B^,  due  to  molecules  which  interact  but  are  free  to  separate 
after  the  interaction.  Equations  are  given  for  determining  each  of  the  three 
parts  of  the  second  virial  coefficients.  A  detailed  treatment  of  these  three 
contributions  is  presented  for  the  square  well,  Sutherland,  and  Lennard- 
JoneB  46  "HSf-potentials.  ^ 


The  mean  lifetimes  of  metastably  bound  double  molecules  are  dis¬ 
cussed  and  it  is  found  that  most  metastably  bound  double  molecules  have 
rn^ati  lifetimes  cunsldersbly  longer  thsn  the  msan  tlm*  between  collisions 
at  ordinary  pressures.  Finally,  an  equation  is  developed  for  the  number  of 
vibrational  levels  of  a  double  molecule. 


* 
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THE  CONTRIBUTION  OF  BOUND,  METASTABLE,  AND  FREE 
MOLECULES  TO  THE  SECOND  VIRLAL  COEFFICIENT 
AND  SOME  PROPERTIES  OF  DOUBLE  MOLECULES 

I.  INTRODUCTION: 

The  equation  of  state  of  a  dilute  gas  can  be  expanded  in  the 
virial  form 

pv/RT  *  i  i-  B<TVv  +  ccryvx  +  ...  (1) 

Here,  B{T),  the  second  virial  coefficient  can  be  expressed  in  terms 
of  the  energy  of  interaction  between  pairs  of  molecules.  If  the  inter- 
molecular  potential,  0  (r),  depends  only  on  the  separation,  r,  and 
not  on  the  relative  orientations  of  the  two  molecules,  then 

B(T)  *  27T N^°[)  (2) 

This  relationship  is  easy  to  derive  from  statistical  mechanical  argu- 

1 

meats  and  has  been  used  since  1908  and  maybe  earlier.  The 
simplicity  of  the  second  virial  coefficient  in  term-  of  the  intermoiecu- 
lar  potential  arises  from  the  fact  that  the  equation  of  state  is  an 
equilibrium  property. 

The  second  virial  coefficient  can  be  broken  up  into  three  parts, 

BfT)  *  a*m  t  Bb(T)  r  BW(T)  (3) 

Here  BJT)  arises  from  collisions  between  free  molecules;  B.  (T)  is 

:  D 

related  to  the  equilibrium  constant  for  the  formation  of  bound  double 
molecules  (or  dimers)  in  the  gas;  and  Bm(T)  is  related  to  the  equili¬ 
brium  constant  for  the  formation  of  metastable  double  molecules  (or 


L.  S.  Ornsteio,  Thesis,  Leiden  1908  , 


1 


2 


dimers)  which,  according  to  quantum  mechanics,  can  dissociate  but 

2 

which,  according  to  classical  mechanics,  are  firmly  bound.  Terrell  Hill 

has  derived  the  equation  for  Qb(T),  the  partition  function  for  bound 

double  molecules,  and  from  this  we  obtain  B,  (T).  In  the  present 

b 

paper,  we  also  derive  the  relations  between  B^T),  B^(T),  and  the 

intermoleeular  potential.  Since  B,  (T)  and  B  (T)  are  related  to  the 

b  m 

equilibrium  constants  for  dimerization,  a  knowledge  of  these  quantities 
is  useful  in  developing  a  theory  for  the  variations  of  the  transport 
coefficients  with  pressure  due  to  the  formation  of  dimers.  The  dissoci¬ 
ation  of  metastable  molecules  and  its  effects  on  the  transport  proper¬ 
ties  of  a  gas  will  be  discussed  Later.  It  will  be  seen  that  the 
dissociation  of  metaetable  molecules  leads  to  a  small  quantum  cor¬ 
rection  for  the  transport  coefficients  at  law  temperatures. 

Whenever  the  potential  is  purely  .repulsive,  so  that  no  bound  or 
metastable  dimers  exist,  the  second  virial  coefficient  can  be  expressed 
in  terms  of  a  collision  integral  which  involves  the  angles  of  deflection 
of  the  binary  collision  trajectories.  Detailed  equations  are  given  for 
determining  Bf(T),  B^tT),  and  D^IT)  for  an  arbitrary  intermolecu- 

lar  potential.  For  molecules  satisfying  the  square  well,  Sutherland 

-6 

(attractive  term  proportional  to  r  ),  or  Lennard- Jones  (6-12) 

t  otential,  tables  of  values  and  graphs  arc  given  for  B*  =  B,/b  , 

t  f  o 

h*  -  B.  /b  ,  and  B*  =  B  /b  as  functions  of  the  reduced  temperature 
b  b  o  m  m  o  r 

>»■  3 

T  =  kT/€  .  Here,  b  =  (2/3)  7T  N O'  The  novelty  of  our  treatment 

o 

consists  in  distinguishing  between  the  bound  and  the  metastably  bound 
double  molecules.  The  introduction  of  the  two  types  of  double  molecules 


T.  L.  Hill,  J.  Chem.  Phys.  617,  (1955);  "Statistical  Mechanics  , 
(McGraw-Hill,  19S6),  Ch.  5. 


3 


is  necessary  to  remove  the  ambiguities  in  the  various  definitions  of 

2 

double  molecules.  Some  authorities,  such  as  Terrell  Hill  ,  have  not 
considered  the  existence  or  effect  of  the  metastable  species. 

Non-equilibrium  properties  of  gases  such  as  the  transport 
coefficients  are  much  more  complicated  functions  of  the  intermolecu- 
iar  potential  than  are  the  second  virial  coefficients.  For  example,  the 
transport  coefficients  in  the  limit  of  low  pressure  can  be  expressied^ 
in  terms  of  the  collision  integrals 

_n5‘>m  -  ir'y-frri"*  C <«, 


Here  /u  is  the  reduced  mass  of  the  colliding  molecules, m*  ; 

Cm,* 

g  is  the  initial  relative  velocity  of  the  two  molecules;  b  is  the  impact 
parameter  or  the  distance  of  closest  approach  of  the  two  molecules  if 
their  relative  trajectory  persists  in  a  straight  line;  and  X  (b,  g)  is 
the  angle  of  deflection  of  the  relative  trajectories  which  results  from 
the  inter  molecular  potential  (see  Fig.  1)  .  The  angle  of  deflection 
can  be  expressed  in  terms  of  the  inter  molecular  potential, 


X(b,9)»7T  -lbQ  rfc 


(5) 


Here  r^  is  the  distance  of  closest  approach  in  the  actual  collision, 

or  r  is  the  largest  value  of  r  for  which 
m 


I  -  -  it  -  0  <* 

If  0(r)  corresponds  to  a  purely  repulsive  potential,  as  shown  in  Fig. 
la,  then  %  (b,  g)  is  always  positive;  if  0(r)  is  a  purely  attractive 


J.  O.  Hirschfelder,  C.  F.  Curtiss,  and  R.  B.  Bird,  "Molecular 
Theory  of  Gases  and  Liquids",  (John  Wiley,  1954),  p.  4S4. 


4 


Fig.  la.  Trajectory  of  Collision  with  Repulsive  Potential. 

4 

Here  %  is  positive. 


Fig.  lb.  Trajectory  of  Collision  with  Attractive  Potential. 
Here  X  i*  negative. 


potential,  as  shown  in  Fig.  1  b,  then  X  (b,  g)  is  always  negative;  and 
if  ^ (r)  corresponds  to  attraction  for  large  separations  and  repulsion 
for  small  separations  (as  is  usually  the  case),  %  (b,  g)  is  positive  for 
•ome  values  of  the  parameters  and  negative  for  others. 

It  is  convenient  to  define  the  effective  potential  energy  as 

4*  L/rx  <7) 

2 

Here  L  =  Kb  where  K  is  the  initial  relative  kinetic  energy  or  the 

2 

total  energy  in  the  center  of  mass  coordinate  system,  K  /Z. 

Insofar  as  the  separation  between  molecules  as  a  function  of  time  is 
concerned,  the  kinematics  are  the  same  as  for  the  one  dimensional 
motion  of  a  particle  of  mass  and  energy  K  moving  in  the  potential 
field  *•»)  .  The  term  L/r  is  known  as  the  centrifugal  poten¬ 

tial.  The  angular  momentum  of  the  collision  system  is  M  =  jll  gb  and 

2 

the  moment  of  inertia  of  the  collision  system  is  I  *  /t.  r  .  Thus  the 

2  2 

centrifugal  potential  L/r  n  M  l?A,  which  is  the  usual  form  for  the 
energy  of  rotational  motion.  Fig.  2  shows  a  set  of  effective  potential 
energy  curves  for  the  Leonard -Jones  {k  12)  potential, 

<t(r )=  *«[  (f)'1  3  <8> 

Here  (-  )  is  the  maximum  energy  of  attra^  tion  between  two  mole¬ 

cules  and  AT  ia  the  low  velocity  collision  diameter,  it  is  seen  from 
Fig.  2  that  $e£f  baa  an  inflec  tion  point  at  .  8  46  when  L  -  Lc  * 

2.  4624€  <T  .  The  value  of  0  ^  at  it#  inflection  point  is  known  as 
the  critical  energy,  1C  .  The  value  of  K  depends  )upon  the  func- 
tional  form  of  the  inter  molecular  potential.  The  critical  value  of  b  is 

defined  as  b  -  (L  /K  .  Thus  for  the  Lennard- Jones  (6-12)  poten- 
c  c  c 

1  /S 

tiatl,  bc  «  (3/5  )  0"  =  1.754410".  For  values  of  L  greater  than 

Lc»  the  effective  potential  curves  have  no  minima  and  vary 


monotonically  with  the  separation.  For  values  of  L  smaller  than  L  , 
the  effective  potentials  have  both  a  minimum  and  a  maximum.  The 
energy  of  the  maximum  varied  between  zero  and  0.8€  ;  and  the 
energy  of  the  minimum  varies  betweeen  {-€  )  end  0.8#  .  Maxima 
in  p  occur  when  the  inter  molecular  potential  is  attractive  at  long 
ranges  and  the  attractive  term  varies  inversely  with  a  power  of  the 
separation  greater  than  two. 

In  Fig.  3,  there  are  two  cross-hatched  areas,  one  labeled  "B" 
and  the  other  "M".  The  region  "R"  corresponds  to  those  two  mole¬ 
cule  systems  where  the  total  energy  is  less  than  the  energy  of  the 
separated  molecules.  The  systems  in  "B"  correspond  to  bound 
molecules  where  the  molecules  can  only  be  freed  by  a  collision  with 
another  molecule.  However,  the  '*M"  systems  have  energy  greater 
than  that  of  the  separated  molecules.  From  the  point  of  view  of 
classical  mechanics,  these  molecules  can  only  be  freed  by  a  collision 
with  another  molecule,  but  from  the  standpoint  of  quantum  mechanics, 
the  "M"  systems  are  metastable  double  molecules  which  can  dis¬ 
sociate  by  leakage  through  the  energy  harrier.  If  the  half-life  for 
dissociation  is  greater  than  the  average  time  between  collisions,  the 
"M"  systems  behave  like  bound  double  molecules.  But,  if  the  time 
between  collisions  is  long  compared  to  the  half-life  for  dissociation, 
the  "MM  systems  behave  more  like  free  molecules.  The  mean  life¬ 
times  of  metastable  molecules  will  be  discussed  in  section  V.  For 
example,  in  the  c«*e  o i  argon  molecules,  it  is  found  that  the  mean 
lifetime  may  vary  between  .about  10t  **  sec.  to  well  over  lO^sec. 
depending  on  the  particular  value  of  the  relative  energy  and  angular 
momentum  of  the  metastable  pair. 

Now  consider  one  of  the  effective  potential  curves  for  a  small 
value  of  L,  as  shown  in  Fig.  3.  This  effective  potential  has  a 
maximum  value  F{L).  Since  F  is  a  monotonically  varying  function 


of  L,  we  can  invert  this  relation  to  obtain  L  as  a  function  of  F. 

Thus  the  family  of  effective  potential  curves  with  F  less  than  or  equal 

to  K  can  be  characterized  by  their  value  of  F  rather  than  by  their 
c 

value  of  L.  As  is  shown  in  Fig.  3,  the  intermolecular  separation  at 

the  hump  or  maximum  in  the  effective  potential  is  r^(F).  Also, 

r^(F)  is  the  smaller  separation  for  which  the  effective  potential  has 

the  value  F.  In  the  limit  that  F  approaches  zero,  r^  approaches 

infinity  and  r  approaches  o'  .  In  the  other  limit  that  F  approaches 
•  « 

K  ,  r.  becomes  equal  to  r  .  The  distance  of  closest  approach,  r  , 
c  h  t  m 

of  two  colliding  molecules  can  be  characterized  by  K  and  F.  Con¬ 
siderable  interest  is  attached  to  those  collisions  in  which’  K  is  very 

nearly  equal  to  F.  In  such  collisions,  the  two  molecules  orbit  around 

4 

each  other  a  number  of  times  before  they  finally  separate  .  If  K  is 
just  slightly  smaller  than  "  then  r^  is  slightly  larger  than  r^; 
wh.’reas  if  K  is  just  slig'utly  larger  than  F,  then  r^  is  slightly 
smaller  than  r^..  Since  small  kinetic  energy,  K  £  K^,  is  required 
for  orbiting,  it  follows  that  such  collisions  are  only  frequent  at  low 
temperatures. 

There  is  still  another  way  that  we  can  characterize  the  effective 

potential  energy  curves  which  have  humps.  Let  us  define  b^(F)  as 

the  value  of  b  for  which  the  initial  kinetic  energy  K  is  equal  to  F. 

2 

Since  L{F)  ~  Kb  ,  it  follows  that 

b/F)  =  [L(  F)/  F]/Z 

Here  bf  varies  monotonies  11  v  with  F  from  the  value  infinity  when 

F  -  0  to  b  when  F  -  K  .  Thus  b,  might  be  used  to  characterize 
c  c  t 

the  effective  potential  curves  with  low  values  of  the  angular 


4 


See  Reference  3,  M.T.G.  L. ,  pages  45  and  553 


***£•  The  distance  of  closest  approach  r  as  a  fane 

m 

tion  of  b  for  a  fixed  value  of  the  initial  kinetic 
energy  K  less  than  K  . 


momentum.  The  significance  of  b^  is  shown  in  Fig.  4.  Here*  for  a 
particular  value  of  K  (less  than  K£)  the  values  of  r^  are  plotted 
as  a  function  of  b.  From  Eq.  (6)  it  follows  that 

bi=e[i-t£>]  <9) 

When  b  =  b,(K)  there  is  a  discontinuity  in  r  .  When  b  is  just 

slightly  greater  than  b^(K),  the  energy  of  the  system  Is  not  great 

enough  to  get  over  the  potential  hump  and  r^  is  slightly  greater  than 

r  (K).  Whereas,  when  b  is  slightly  less  than  b,(K),  then  r  is 
h  i  m 

slightly  less  than  r^(K).  When  b  is  aero*  r^  becomes  equal  to 
r^(K)  which  is  defined  in  terms  of  the  equation, 

(p(r.)  *  K 

The  meaning  of  r^  and  b^  will  be  explained  in  a  later  section.  How¬ 
ever,  since  0  (r  )  becomes  small  as  r  becomes  large,  it  follows 

m  m 

from  Eq.  (9 )  that  the  value  of  b  approaches  the  value  of  r  as  r 

m  m 

becomes  large  for  a  fixed  value  of  K. 

The  Second  virial  coefficient,  B(T),  for  spherically  symmetric 
molecule's  can  be  written  as 

B(T)  *  -NA*  Qx/V  ♦■NA*Q)VaV  <10} 

where  Qj  and  are  partition  functions  for  one  and  two  particles 

respectively,  N  Is  Avogadro's  number,  V  is  the  volume,  and 
2  2 

A.  =  h  /2TT  m  k  T.  The  partition  function  involves  integrations 
over  the  volume  in  which  the  particle  is  contained  and  over  the  momen¬ 
ta  of  the  center  of  masa  of  the  particle.  The  two  particle  partition 
function  ,  involves  an  integration  over  the  volume  in  which  the 

two  particle  system  is  contained,  an  integration  over  the  relative 
coordinates  of  the  two  particles,  an  integration  over  the  momenta 


of  the  center  of  meet,  end  an  integration  over  the  relative  momenta  of 
the  two  particles. 

Systems  "BM  and  “Id”  of  Fig,  3  have  partition  functions 
Q  and  Q.  respectively.  They  involve  integrations  over  the  same 

variables  as  Q^.  However,  the  integration  limits  are  restricted  so 
that  the  two  molecules  remain  in  the  regions  *9*'  and  *,M"  indicated 
in  Fig.  3.  The  integration  limits  are  discussed  in  more  detail  in 
Section  III.  Then,  B^(T)  and  BJT)  are  deflaed  by  the  equations 

* -n  a4q„/v  mi 

? 

B»(T)«-U*Qw/V  <u> 

The  partition  function  of  all  other  possible  systems  in  Fig.  3  is  Q2f. 
Again,  contains  integrations  over  the  same  variables  as  Q ^  , 

but  the  integration  limits  are  such  that  the  two  molecules  remain  out- 
side  of  the  regions  "B"  and  Mht" ,  One  sees  that  these  systems  cor¬ 
respond  to  molecules  which  interact  but  are  free  to  separate  to 
infinite  distances  after  the  interaction.  The  definition  of  B^(T)  is 

^  (T)  «  -  N  A*  ■»  NA‘  Q  *  /i  V  (13) 

The  role  of  double  molecules  in  the  second  virial  coefficient  has  long 
been  recognised.  We  will  now  show  that  the  equilibrium  constant 

tt(T)»  <M) 

1.  fim.  ky 


am  -  *  s»m 


(15) 


In  Eq.  (14),  and  n^  arc  the  number  of  moles  of  single  and  double 
molecules  (  i.e.  n^  =  n^  +  n^m  ;  the  subscripts  b  and  m  refer  to 


bound  and  metastable  double  molecules  respectively)  and 


t  j  and 


are  activity  coefficients.  At  low  pressures,  the  activity  coef- 

2  5 

ficients  approach  unity.  From  the  equations  given  in  Hill's  paper  , 

one  easily  obtains  the  relation 


*tV/n*  *  -A?  (,Qlb«-Qx~)N/V  4  £71  *♦*•••  (16) 

The  right  side  of  this  equation  is  a  power  series  in  n  which  is  defined 
by 


n  *  n,  t  ♦  an h 

The  number  of  moles  of  "triple  molecules"  is  n^f  etc.  The  coeffi¬ 
cient  "a"  involves  integrals  over  certain  regions  of  the  phase  space 
of  one,  two,  and  three  molecules  respectively.  Higher  coefficients 
depend  on  integrals  over  part  of  the  phase  spice  of  greater  numbers 
of  molecules.  Comparison  of  Eq.  (16)  with  Eqs.  (11)  and  (12) 
shows 


46^)  -rant-...  i1 

2 

The  equilibrium  constant,  X(T),  is  equal  to  n^V/n^  in  the  limit 
that  n  approaches  zero  because  under  this  condition  the  activity 
coefficients  approach  unity.  Thus,  Eq.  (15)  forJ^(T)  follows  from 
Eqs.  (14)  and  (17). 

The  equilibrium  constant,  %  (T),  can  be  written  as 


5 


Reference  2. 


(18) 


xm«xkm  +  xmcn 

wb«rt  ^^(T)  is  ths  equilibrium  constant  for  bound  do  obi*  moUcuUs 

sod  *  (T)  is  the  equilibrium  constant  for  me  testable  double 

m  * 

moUcuUs.  It  follows,  for  reasons  similar  to  thoso  givoo  tor  K  (T), 
that  -Xfc<T)«B.<T)  »■*  -X„(T)«B i  (T)  . 

D  D  mm 

In  section  IH  it  will  be  seen  that  B.  (T)  and  B  (T)  must  be 

d  m 

negative  so  that  the  equilibrium  constants  are  positive  as  they  should 

be.  This  is  to  be  contrasted  with  the  results  of  previous  attempts  to 

relate  the  second  virial  coefficient  to  the  equilibrium  constant  for 

dimer  formation.  In  these  previous  attempts,  because  oi  the  approxi- 

6 

mat  ions  made,  either  a  negative  equilibrium  constant  is  obtained  at 

moderate  and  high  temperatures  or  the  introduction  of  an  excluded 

volume  due  to  the  finite  else  of  molecules  obviates  the  necessity  ri 

7  8  9 

dealing  with  a  negative  equilibrium  constant  *  '  .  The  equations  de> 
rived  in  section  Ill  will  relate  the  equilibrium  constant  to  the  inter  - 
molecular  potential  without  the  introduction  of  any  empirical  concepts. 

The  effects  of  interactions  between  two  free  molecules  is  repre¬ 
sented  by  Bg(T);  Eq.  (55)  shows  how  Bf(T)  depends  on  the  angle  of 
deflection  due  to  s  collision.  If  it  is  desired.  B^(T)  can  be  related  to 
the  empirical  concept  of  excluded  volume  by  supposing  that  a  gas  is 

*  H.  W.  Woolley,  J.  Chem.  Phys,  21,  236  (1953). 

7  X.  O.  Hirschfilder,  T.  T.  M&lure,  and  I.  F.  Weeks,  J.  Chem. 
Phye.  ifl.  201  (1942). 

8  W.  We  liner,  Jr.,  J.  Chem.  Phye.  ££,  153  (1954). 

*  R*  Ginell,  J.  Cham.  Phye.  2395  (1955). 
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composed  of  n^  mole*  of  single  molecules  and  moles  of  double 
molecules  satisfying  the  equation  of  state 

P(V~k>)  *  (71,  i-tiJKT 

Except  for  the  excluded  volume,  b',  this  corresponds  to  the  perfect 
gas  equation.  The  number  of  single  and  double  molecules  are  deter¬ 
mined  by  Eq.  (14)  with  the  activity  coefficients  set  equal  to  unity.  If 
we  take  V  to  be  the  molal  volume,  (regarding  the  substance  to  be 
composed  of  monomers)  then 

71,  +  ZVX  *  I 

because  the  concept  of  molal  volume  is  based  on  the  notion  that  all  the 
molecules  are  Single.  From  the  last  equation  and  Eq.  (14),  it  follows 
that 

For  large  volumes,  expanding  the  right  side  of  this  equation  in  powers 
of  1  /  V  gives 

TV  =  X(T)/V  -  HfXCOAV)*' 

Thus 

PV/RT  -  [J  -  nj  (I  -  b'/vr' 

=  [I  -  X/V  +H(X/vf  -...111  +  b'/V  +  ( b7v)‘  + ...] 

=  I  +  Ib'-KW  +  (b'*-b'X  +HX‘yv‘+... 


so  that  the  second  virial  coefficient  becomes 


B(T)  -  b'  -  X(T) 

The  excluded  volume,  b',  can  then  be  identified  with  B,  (T). 
Ordinarily,  B^(f}  in  positive  corresponding  to  positive  values  of 
the  excluded  volume*  However,  at  low  temperatdres  Bf(T  ),  sad 
therefore  Ole  excluded  volume,  becomes  negative!  Thus,  this 
simplified  analysis  shows  the  limitations  of  the  concept  of  excluded 
volume. 

We  will  now  proceed  with  the  derivations  of  the  relations  for 
B^(T),  B^(T),  and  Bm(T).  Actually  these  relations  were  obtained 

when  we  tried  to  extend,  to  molecules  obeying  an  arbitrary  inter- 
molecular  potential,  a  result  which  R.  Byron  Bird10  has  obtained 
for  molecules  obeying  a  monotonically  varying  inter  molecular  poten¬ 
tial.  Bird  had  us#<jf  quantum  mechanical  arguments  to  express  the 
second  virial  coefficient  in  terms  of  the  Sag^s  of  deflection  of  the 
trajectories.  For  monotonically  varying  intermolecular  potentials, 
no  bound  or  mOtaetable  double  molecules  can  exist.  Our  results, 
obtained  by  classical  mechanical  arguments,  agree  with  Bird's  for 
this  special  type  of  potential. 


II.  BIRD'S  DERIVATION  OF  THE  SECOND  VIRLAL  COEFFICIENT 


FOR  MONOTONICALLY  VARYING  INTERMOLECULAR 
POTENTIALS10: 

The  quantum  mechanical  expression  for  the  second  virial  coef¬ 
ficient  of  a  gas  of  spherically  symmetric  particles  in  the  case  where 
the  potential  is  monotonically  decreasing  with  increasing  separation 
of  the  particles  is** 

BCf)  *  i  -  2mNA’  £(l-t  +l)(V* <  w> 

2,  .  2 

where  (X)  is  the  phase  shift  and  X  =  2/VK/Jt  .  In  the  case 
of  Fermi-Dirac  statistics,  the  plus  sign  is  used  aftd  the  sum  over  Z 
includes  only  odd  values  of  A  .  For  Bose -Einstein  statistics,  the 
minus  sign  appliea*  and  the  summation  is  over  even  values  of  A  .  In 
the  case  of  Bpltamann  statistics,  B(T)  is  the  average  of  the  Fermi  - 
Pirac  and  Bose-Einstein  results.  Thus, 

«T)  -  - 2-NA*  5  (U  ,20, 

where  the  summation  now  extends  over  all  values  of  A  . 

To  obtain  the  classical  mechanical  expression  which  corre¬ 
sponds  to  Eq.  (20),  the  following  steps  are  carried  out:  (1)  The 
integral  in  Eq.  (20)  is  integrated  by  parts.  Use  is  made  of  the  fact 
that  7(^(0)  s  0.  (2)  The  sum  over  A  is  replaced  by  an  integra¬ 

tion  over  JL  .  The  result!*  - 

m  •  •  • 

>0 

Reference  3,  M.  T.G.  L. ,  p.  51.  We  wish  to  thank  R.B.  Bird  for 
making  available  his  unpublished  notes  giving  his  derivation. 


11 


Reference  3,  M.  T.G.L. ,  p.  409 
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sen  »  -  {2i. 

Let 

A'lJU+i)!*  =  AS^  (22) 

&K  *  /*$  (23) 

Through  the  use  of  the  W.K.B.  method,  it  is  known  that  to  the  first 
approximation1 2, 

*<b,g)  ”  k  sS‘  1241 

Eq.  (24)  becomes  correct  in  t1  e  limit  of  large  values  of  K  or  / 
where  the  correspondence  principle  requires  quantum  and  classical 
mechanical  results  to  be  the  same.  It  should  be  noted  that  Eq.  (24)  is 
only  valid  if  there  is  only  one  turning  point,  and  the  inter  molecular 
potentiates  monotonies  11  y  decreasing  with  increasing  separation  of 
the  particles.  Integration  of  Eq.  (24)  gives 

\  *  j[*  (25) 

W  nsui  it*.*  { *  J  r  i)<ii  >■  A  A*-  ^  but,  tw  utittti,  muii  u«  Uv»’  t»i 

integration  over  b  and  b1  is  changed,  Eq.  (21)  becomes 

am-  2N7r^rAT)’,Af e'^^tfZh'JhJdK  m 


12 


Reference  3,  M.  T.G. 


L. ,  Section  10.3b,  p.  687. 


This  is  the  equation  Bird  obtained  and  is  valid  in  classical  mechanics 
if  the  potential  is  monotonic  decreasing  with  increasing  sepatation  of 
the  interacting  particles.  Below  Eq.  (45),  it  will  be  shown  by  f.tnrti 
classical  mechanical  arguments  that  when  the  potential  is  monotonic. 
Eq.  (26)  is  an  alternate  means  of  writing  the  more  usual  expression 
for  the  second  virial  coefficient,  i.  e.  Eq.  (2). 
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IU.  DERIVATION  OF  THE  SECOND  VIRLAE  COEFFICIENT  FOR 
ARBITRARY  INTERMQLECULAR  POTENTIALS  IN.. TERMS  OF 
THE  DOUBLE  MOLECULES  AND  COLLISION  PROCESSES: 

In  order  to  investigate  the  contributions  to  the  second  virial 
coefficient  from  the  bound  and  metastably  bound  double  molecules,  it 
is  of  interest  to  relate 

B'rn  -  2NTr,/VATfnT l27) 

O 


to  B(T)  as  given  by  Eq.  (2)  for  the  case  where  the  potential  is  not 
necessarily  monotonic.  Only  in  the  case  of  the  monotonic  potential  is 
B'(T)=B(T). 


>1  J  >tb^  db  ,  is  found  by  taking  the  limit  of 


The  intej'.r  ■ 

2 

J(b  db  as  b  approaches  infinity.  Let  us  define  r  an  the  value 

A  1 


of  r  corresponding  to  b  and  the  given  value  of  K.  If  b  is 
ml  1 

sufficiently  large,  then  r^  is  very  nearly  equal  to  b^  irrespective 
of  the  value  of  K.  Using  Eq.  (5)  and  the  definition 


2  -2 


2,^1  ~h 


G(b.  r)  =  b~  r  [l  =  0{r)/K-bt*/r*“J  we  ootam 

f^b'dlb  =  TTbjVs  ”  )  G(b,r)drdb% 

If  the  integration  over  r  is  performed  last,  Eq.  (23)  becomes 
_ ,%JL  rr*(K)  .. 


(28) 


J^Xb'ib  -  Tbf/3  -/Vl"'fltr  /'  ''  *  JG(b fr)db" 

~  /‘‘"’dr  Glb/lab1  -  f' 

—  otb  f* G(b,r)dr  (2< 

*  r. 
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When  use  is  made  of  the  fact  that 


jGCb/Ub*'  =  -2[-§rt('“  ~  j^)  +  bX(l-  ^p  ”  ] 


and  some  of  the  resulting  terms  are  combined  into  one  integral, 
Eq,  (29)  can  be  written  as 

=  1Tb’/3  -  $fr‘  r'[i-*f?]Adr 

•  i;(K) 

+  2K"'/iC(K)  -  )o'  *£  j°°G(b,r)c(.r 

where 


C(K)  = 

fk(K,[-|-r1(K-#m-^)“  -  K^(K-«fr)-  0  f]d.r  (32) 

ff(K) 

The  last  integral  in  Eq.  (31)  will  be  evaluated  in  the  limit  that 

*  1  .  ~ 

Tj  (or  bj) — >  oo  .  Let 


J  =  (b'dbl  f  G(b, rUr 

0  r. 


,b«r* 

rbt 

,6'*f 

=  [  '  j 

+  1 

+  J 

b  *  0  V, 

b  *  r' 

k'*r,x  bi 

I 
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nr 

First,  an  upper  limit  for  J  will  be  calculated,  and  this  will  be  shown 
to  He  equal  to  the  lower  limit  for  J  as  r^  — »oc  , 


If  the  positive  quantity  -  0(r)/K  is  neglected,  an  upper  limit  for  I 
is  obtained.  Since 


jV*o-£rv 


(I 


i<rr'(h/0  f 


I„,  -  -  2 r* ( t/6  -  2*)  (34) 

Similarly, 

a*.,  -  2  jVsin-,(b/b1m  - 
2b^[w  sm-ft/u)  -  o/»tX('c/b.f+2}0-^/b,)‘}v*J  <35> 

Use  of  Eq.  (S0j  gives 


If  the  negative  terms  of  Eq.  '  3o)  are  neglected  and  -  ^(r)/K  is 
replaced  by  tbs  larger  quantity,  ■  $Xt  j)  K,  one  finds  that 


(36) 


2/vrto-^r 


-  *0 -tyfti-rn  t  a^fi-^Prfh-O 

One  finds  ahd  Ilmax  approach  aero  as  r  ^  — *►  oo  . 

case  of  11^^,  the  assumption  must  be  made  that 


(37) 


In  the 
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lim  r*  *o  (38) 

r-^o* 


Thus, 


lim  Jwax  =  TTr;V3  “  vr.Vs  (39) 

A  lower  limit  for  J  is  found  by  replacing  -  0(r)/K  in  Eq.  (33)  by 
the  larger  quantity  -  $( r  )/K. 

if  b^fcfbr- -  zr*[, -%>]*[$  -  (40, 

lim  J-,.  «  THjVs  -  Mr; •/<♦  (4i, 

3  3 

Eqs,  (39)  and  (41)  show  that  lim  J  =  TTr^  /3  -  4r^  / 9.  Eq.  (31) 

r  —•ao 

1 

becomes  (on  taking  the  limit  b  ►  oo  ) 


lim  j  1  Xb  db  —  lim 


TTb’/s  -TTCV3  1*  9^/9 


-$/'  ]"*<ir  tiK'^crw 

fyw 


=  to*  { 

b,-«. 


-  I  fr'  rio-^pf*  +  iK'**C(K) 

*  Vkj  l 


(42) 


Thus 


fzkVI,  -  +aR— cw 


In  both  Eqs.  (42)  and  (43)  use  has  been  made  of 


(43) 


i 


=  *\  (1  -  ^(r  )/K)*  and  lim  r  =  b  . 
1  1  1  fc^-aao  1  1 


>■ 


Therefore, 


e'rn  =  e  **V*{  +2^00}  (44, 

It  ia  desirable  to  perform  the  integration  over  r  in  T^q.  (44)  after  the 
integration  over  K.  Fig.  5  shows  the  region  of  integration.  The 
distance  of  closest  approach  when  ^eff  =  0 (r)  (L  =  0  in  Eq.  (7))  is 
O'  when  K  =  0. 


In  the  case  of  monotonic  decreasing  potentials,  C(K)  is  aero 

because  r.(K)  *  r  (K)  for  ill  values  of  K.  Also,  in  this  case, 
t  n 

Eq.  (45)  would  be 


because  <JT  *  0©  .  After  some  manipulation,  it  is  found  that 

£*£V(rK-*r>}M  -Kv‘)e*K'ATdK<ir 

*  f  ff w? e'**7 dK 

so  fhfct  Eq.  (4*1)  is  Identical  yrith  Eq.  (2).  Thus,  the  validity  of  Eq. 
(26)  ia  classical  mechanics  is  demonstrated  without  the  use  of 
quantum  mechanics. 

When  Eq.  (45)  is  used,  Eq.  (44)  becomes 
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B-m  - 


°  ® 

t  lfQ‘e'MTK*tdKdr  -  $) -K**Jdl Kir 


(Kr) 

rK* 


[♦  2  l  tC(K)e‘**TdK 


(46) 


In  the  last  term  of  Eq.  (46),  the  upper  limit  for  the  integration  over 

K  **  *  becau*e  for  values  of  K  above  K  ,  r  (K)  =  r  (K). 
c  c  l  h 

Then  ,  by  using 

£jK-pto]**c','*rM  »  -f  ir'*(AT)t/1e*(r>9(T 

ffrV^VW  =  i-T*aTf*frdr  -  l°e'KKTK**[toof dK 

•01  f  h  o 

£VV**T4K  =  $  TT  Uk(ATf/’- 


and 


(VftrjrV^K  =  ln*(JCrFe«rVAT -(ATre^nir®) 


one  finds 


b'(T)  «  iNlT^AT) 


M 


[  +  i^COOf^dK 

^*r*  I  (n»  *)  i*  »n  incomplete  gamnu  function. 


V^(AT n~r'[i-  e-*WATJdr 

t  $  UiT)’/1fr'e*'''*rr(ij-$p)dr 


(47) 


It  is  known  that  the  partition  function  for  bound  pairs  of  mole- 
.  .  13 

cules  xs 

Qlk  =  znVA-'iSeMTm,-$ym<ir 

Thus,  Eq.  ( 1  1)  gives 

B  (T)  =  -2N7T  (})]dr  (48, 

o  (T 

The  use  of  Eqs.  (2),  (47),  and  (48)  together  with  the  equation 

m,-|^)/r(S)  =  r(lr%p)/rti)  -  e*mT[ -$$f/r(i) 

gives 

B(T)  =  E^IT)  t  B'(T)  +  •f‘N7T'/it  fl[’  *r 

-  jf*  C(K)  e-H/*T  *LK  (49, 

In  order  to  find  the  contribution  of  the  metastably  bound  double 
molecules  to  the  second  virial  coefficient,  we  start  with  the  equation 
for  the  partition  function  of  two  particles: 

*  2^*le  /*'Tdxdy<iz<l$d0d.rdpK  dp,  dp9  4p^  dpr 

H  *  rpf  hpJ)/Mm  i-  (pf  t  p£/si*xe)/v\rx  +  prVw  i  <p(n 

Here,  x,  y,  and  z  refer  to  the  cartesian  coordinates  of  the  center 

of  mass  while  the  other  coordinates  are  the  relative  coordinates  of  the 

2  - 

two  molecules.  If  the  change  of  variables,  P#  =  p#  /(mr  kT.  2  , 

2  -  4 

=  p^  /(mr  kT)*,  and  =  p^/fmkT)*  is  made,  and  the  inte¬ 
grations  over  x,  y,  z,  B  ,  $  ,  p  ,  p  ,  and  p  are  carried  out/ 

XV  z 
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Reference  2. 
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Q.  -  e«p[-  XT  *]drdP'dPf  dpr  (50) 

The  total  energy  in  the  center  of  mass  system  is  given  by 

K  -  AT  (P*  i -Pf)  +  4>  (r)  +  ATP,.1' 

Mao, 

(rx  +  P,*)AT  -  Kb Vrx 

2  2 

where  Kb  /r  is  the  centrifugal  potential.  Only  the  sum  of  the 
squares  of  and  is  important  in  determining  the  centrifugal 

potential,  so  a  change  to  polar  coordinates  is  made  using 

ff  dP0  dP4  *  T  (AT)'1  jV'Vf'Kb*) 

When  the  variable  is  replaced  by  the  variable  K, 

*  VlT^CATr^X* H4kTal(Kbx)dtdK  (51) 

To  obtain  B^T),  the  integration  in  Eq.  «'51)  must  be  restricted  to 

the  regions  indicated  in  Figs.  6a  and  feb.  For  a  fixed  K,  Fig.  6a 

? 

shows  the  effective  potentials  for  the  two  extreme  values  of  Kb 
which  are  allowed.  The  lower  of  these  is  Kb^  and  the  upper  is  Kb^. 
Fig.  shows  the  region  of  integration  for  a  fixed  K.  The 
condition 

Kb*  «  rx[K  (22) 

is  true  at  tha  turning  points  Of  motion.  From  its  definition,  Eq.  (12), 


<p  =<t>(')i-Kbl/r*- 

'  C++  T 


Fig.  6a.  The  effective  potentials  when  b  =  and  b  =  b 


Fig.  6b.  The  region  of  integration  indicated  in  Eq.  (53) 
for  a  fixed  value  of  K. 
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VT)  '  "S^P*  i  {,»  L*  fK-^r,-T*t  ]  (S3) 


•o  that 


vr»  -  -  rcr^-^  r  <-» 

Combining  Eq*.  (49),  (54),  and  (3)  show*  that 

qro  -  but)  +  4-  r‘f'^w3vi^r 

-  I"* co°  «"’<MTetK  -«w")  <«> 

Thua,  the  second  virial  coefficient  has  been  divided  into  three 

parts.  Bfe(T)  is  given  by  Eq.  (48),  Bm^T)  bY  E9*  (54),  and  Bf(T) 

by  Eq.  (55).  In  the  case  where  the  potential  is  monotonic  decreasing, 

all  the  terms  except  B’(T)  are  zero  and  B(T)  is  given  by  Eq.  (26). 

For  arbitrary  potentials,  B.  (T)  and  B  (T)  have  to  be  considered 

d  m 

in  addition  to  B^(T).  In  B^(T),  the  terms  added  to  B'(T)  are  a 
result  of  the  fact  that  there  is  more  than  one  turning  point  (see 
comment  below  Eq.  (24)  ). 
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IV.  EVALUATION  OF  THE  THREE  PARTS  OF  THE  SECOND  VIRIAL 
COEFFICIENT  FOR  SEVERAL  INTERMOLEC  ULAR  POTENTIALS: 

The  results  of  the  previous  section  will  now  be  applied  to  the 
square  well,  Sutherland,  and  Lennard- Jones  potentials.  For  these 
potentials,  it  is  convenient  to  write  the  equations  in  terms  of  reduced 
variables.  This  also  permits  an  easy  tabulation  of  numerical  results. 

30C  jjg 

Let  r  =  r /<r  ,  b  =  b/c r  ,  K  =  K/€  ,  <P  =  #/€  ,  T  =  kT/<£  ,  and 
D(K*)  =  (ere**)'5  C(K).  Thus, 


Bf(r-)  *  -3 /W*VT*  USfcrsrJl  oLr* 

&~(T')  =  (*'  M(K")  e-'s'/TbK' 

m  '0 

M(K")  *  00  ,»*[**- dr" 
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(a)  The  Square  Well  Potential. 
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(57) 

(58) 


(59) 


(60) 


This  potential  is  given  by 


.  *  * 

$  (r  )  = 


* 

r  C  1 


* 

<  r 

* 

r 


<.  R 
>  R 


0 
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*  * 

The  expression  for  B,(T  )  is  given  by  Eq.  (56), 

D 

8*rr*)-  a-Tr'^i-*1)  r(Va,  t*"*)  eT*'*  (6i) 

Because  of  the  discontinuity  in  the  potential  at  x  =  R,  some 
special  features  arise.  In  the  integral  for  B*  (T  ),  two  cases  occur 

t- 

depending  on  the  value  of  K  .  The  separation  of  the  two  cases  occurs 
when 


<C<» 


K-bf'/H1 


or 


-  k%;vrx 

♦2  2 

Since  »  R  ,  it  follows  that 


K*  =  (  R  -  / ) 


-i 


In  case  I,  K*  <  (R2  -  l)"1  and  r*(K*}  *  1.  In  case  II,  K*»  (R*  -  if1 
*  *  1 

and  r.(K  )  >1.  Figs.  7a  and  7b  show  the  two  cases.  Also, 

*  * 
because  of  the  discontinuity  of  the  potential  at  r  *  R,  it  follows  that 

* 

K  is  infinite  ini  tend  of  finite  as  is  the  case  for  smoothly  varying 

^  a  a 

potentials.  To  find  B*  (T  ),  the  integration  limits  for  M(K  )  must 

m 

first  be  determined  from  Eq.  (52).  For  case  I,  one  obtains 
r*(K*)  -  1  and  rh*(K*)  -  R.  For  case  II,  rf*(K*)  *  R[K*/(K*  ♦  1)]  l 

*  4-  * 

and  r  (K  )  *  R.  The  integration  over  r  gives 

A 


The  fir  at  integral  in  Eq.  (62)  is  easily  integrated  in  terms  of  the 
exponential  integral 


Fig.  7a.  <p  effversus  r  for  the  Square  Well  Potential. 
Case  I:  K  <  (R2  -  l)-1  and  r*  (K*)  =  1. 


* 


<t>  eff  versus  r  for  the  Square  Weil  Potential. 
Case  II.  K  >  (R  -  1}  and  r *  (K*)  >  1. 
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EA-t) 


(63) 


The  second  integral  in  Eq.  (62)  tan  be  evaluated  by  making  a  change 

2  *-  Z 

in  variables  to  y  -  [  1  (R  •  1)K  j/R  and  then  expanding  the 
exponential  term.  Thus 


where 

B  j*  X%/x(i -x)K'*dx  (65) 

•/*v  ’  ;0 

is  an  incomplete  beta  function. 

The  equations  for  'X.  ,  to  be  used  in  Eq.  (59)  are  gi*'er.  by 
14 

Holleran  and  Hulburt  The  integrations  in  Eq.  (59)  can  be  carried 

out  to  give 

B*  (T*)  *  I  -  ( A*  -  Of  eT*  -  ()  -  bVt")  -  O'*)  (66) 

As  a  check  on  the  work,  Eq.  '3">  'in  its  reduced  form)  can  be  used 
in  Eq.  '66,.  When  this  is  don»  one  obtains  the  well  known  result  for 
the  second  virial  coefficient  of  particles  interacting  with  a  square 
well  potential. 

Table  I  gives  the  valut  a  of  the  coefficients 
which  appear  in  Eq.  (64  .  The  various  contributions  to  the  second 
virial  coefficient  are  shown  in  Table  II  for  three  values  of  the 
parameter  R.  Fig.  8  is  a  graph  for  the  case  R  *•  1.5. 


14 


E.  M.  Holleran  and  H.  14.  Hulburt,  J.  Chem.  Phys.  232  (1951k 
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Table  1. 


The  Square 


Well  Potential  Coefficients 

¥htch 


e  * 

Occur  in  |q.  (64)  for  B  (T  ) . 

m 
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X 
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X 
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io-6 

i 

8 

1.50388 
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(b)  The  Sutherland  Potential. 


The  potential  for  this  model  is 


4>*(r*)  =  oo 

* 

r  <  1 

♦*(r*)  =  -r*"* 

* 

r  >  1 

The  equation  for  B*(T  )  will  be  given  as  a  power  series  in 

*  „  i  D  x 

T  .  An  expansion  for  e  P  (3/2,  x)/T  (3/2)  is  first  developed;  this 
expansion  will  also  be  useful  in  finding  B*  (T  )  for  the  Leonard- Jones 


potential.  The  incomplete  gamma  function  is  a  particular  case  of  the 

15 

confluent  hypergeometric  function  which  is  defined  by 


$ra,c;x) 


f  Cl la  ** 

(c\  nj 


(67) 


where 

<a\  *  i 

(a)  -  a(a+  o...  (a  - o 

Thus, 

x  f  i'<nri)Ul  (68) 

1 r7  ho 

Here,  the  second  step  follows  from  Kummer's  transformation  and  the 
third  step  from  Eq.  (67). 


A.  Erdelyi,"Higher  Transcendental  Functions'1,  Vol.  I, Bateman  Manu¬ 
script  Project,  (McGraw-Hill,  1953).Ch,  VI  of  this  book  has  a  summary 
of  all  the  properties  of  the  confluent  hyper  geo  me  trie  function  used  in 
this  paper.  It  also  has  information  on  the  Whittaker  functions  which 
are  introduced  in  the  discussion  of  the  Lennard-Jones  potential. 


If  the  expansion,  Eq.  (68),  is  used  in  Eq.  (56),  integrals  of  the 
form  f  r-1  ***  ^  <tr*  *  — - - -  result;  Eq.  (56)  becomes 


7**0 


Hn  (  71  +  l )  1 _ 


T»n 

f 


(69) 


.  * 

Eq.  (69)  is  valid  if  Y  >  2.  It  will  be  shown  below  that  B  (T  )  also 

m 

is  defined  for  Y  *  Z.  This  is  to  be  contrasted  with  the  result  for  the 

second  virial  coefficient  for  this  potential  which  is  defined  only  for 

Y  >  3.  The  fact  that  B*  and  B*  are  defined  when  y  >  2  implies 

b  m 

that  an  equilibrium  constant  for  the  formation  of  double  molecules  can 
be  calculated  (see  Eq.  (15))  for  some  cases  where  the  second  virial 
coefficient  does  not  exist. 

For  the  special  case  where  Y  ~  6, 


B*  (T*) 


-  H 


77-  "4 


(70) 


where 


a..  *  ■  2.  JELl .  (7i) 

We  find  that  aQ  «  1.  60666  x  10'\  a  j  =  3.  33333  x  10 “2, 

a,  =  6.  34921  x  10*3,  d,  -  1.  05820  x  10“\  a  =  1.  53920  x  10*4, 

2  3  4 

a_  =  1.97334  x  10'5,  a,  =  2.  25524  x  10"6,  and  a,  =  2.  32157  x  10'7. 
t)  0  7 

*  *  *  *2  * 

To  find  B  (T  ),  one  must  first  determine  K  b,  ,  K  ,  and 
m  f  c 

A 

integration  limits  for  M(K  )»  From  the  form  of  the  potential,  it  is 

obvious  that  rf  (K  )  =  1.  One  finds  K*  ^and  (K*)  « 

h~jj^  A  by  solving  the  equations  d  0#££ /dr  =  0  and 
^f£(r  )  *  K  .  By  letting  r£  (K  )  *  r^fK  )»  *•  found  that  K*s 


are  used  and  the 


*  * Z 

Whe  n  the  above  equations  for  K  b^.  ,  r ^  , 


substitution  t  = 


Vjf 


*  * 

and  r 

n 


is  made,  the  integral  M(K  )  becomes 


”<**)*  («  tr"+|J dt 


(72) 


If  IT  >  2  in  Eq.  (72),  it  is  easy  to  see  that  Eq.  (57)  gives  a  finite 

value  for  B^.  It  should  be  noticed  that  the  polynomial  under  the 

square  root  sign  always  has  a  double  root  at  t  -  1. 

The  integral  of  Eq.  (72)  can  be  evaluated  for  several  choices  of 

Jf  but  unfortunately  a  different  method  must  be  used  for  each  choice. 

2 

We  shall  consider  only  the  case  f  =  6.  When  the  substitution  z  =  t 

2 

is  made  in  Eq.  (72)  and  ( 1  -  z)  is  factored  from  under  the  square 
root  sign,  Jd(K  )  is  easily  evaluated.  Thus, 


M(K*)  -  f'-tltW'H'1  + 


-  coth’1  i"*  +  coth-1  [l  f-  (73) 

On  expanding  the  exponential  in  Eq.  (57)  and  changing  the  integration 

^  1/3  1 

variable  to  y  =  [l  +  (4K  )  '  ]  *  ,  Eq.  (57)  can  be  integrated  to  give 


b*  m 


J  »  9/t. 


(74) 


where 


s'"  (Z‘nt3X*nt3)7i!  <• 


^(sni-q 


r*o  (3 TJ'fl-r)]  r l  (4 Tlf-7-i.r)  J 


(75) 
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We  find  cQ  =  1.  23065  x  10”  ^=-3.  75800x10"?  c^  =  1.  20 161  x  10  2, 

c,  =  -3.  51826  x  10'\  c  .  =  9.  23385  x  10"4,  c  =  -2.  1740  x  10'4, 

3  4  5 

_5  -6 

c  ,  =  4.6196  x  10  ,  and  c_  =  -8.9247  x  10 

6  • 

The  total  reduced  second  virial  coefficient  for  this  potential 

,  . 

when  f  *  o  i» 

B*  (Tn)  -  I  -  £  fin  <76> 


where 


d 


K 


_J _ 

(xn-i)ni 


(77) 


Here,  d^  *  1,  d2  =  1.  66666 x  10‘  \  =  3.  33333 x  10'2,  d4=  5.  93238 xlO"? 

d5  =  9.  25926  x  10"4,  dfe  =.  1.  26263  x  10 "4,  d?  =  1.  52625  x  10 "5,  and 

dQ  =  1.  65344  x  10“6. 
o 


When  Eq.  (3)  is  written  xn  reduced  form,  B^(T  )  can  be  found 

j|{ 

by  subtracting  B*(T  )  and  B*  (T  )  from  B*(T  ).  Table  III  and 

b  m 

Fig.  9  show  the  contributions  to  the  second  virial  coefficient  for  the 

♦  -6 

Sutherland  potential  with  the  attractive  term  equal  to  -r 
(c)  The  Lennard-Jones  (o-12)  Potential. 

The  equation  for  this  potential  is 


p*(r*)  - 


(78) 


16 


Reference  3,  M.  T.G.L. ,  p.  58.  This  follows  from  the  equation 
given  for  B(T)  which  is  valid  for  V  >  3  after  the  typographical 
errors  are  corrected. 
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^  *  ♦  %  * 

To  obtain  Bfa(T  ),  lot  x  =  *  $  (r  )/T  and  ooo  Eq.  (68).  In  tho 
•erics  for  B*  (T  ),  integrals  of  the  form 

D 

-  fr** 

I 

-  i  (f.fl  (  •’*-“  [I  -  r—  f  “*•  d(r‘-*) 

appear.  Rcpoated  integration  by  parts  gives 


a  «  f  X\tf*  a**'  nl 

'  It*/  6  ( 2  n  j*X2n  ♦'iXan  *91. .  (h  n 


Thus 


fif ct*)  *  7T,^'t**/*  ^  f** 


whore 


^  (*jiTs)l -  vou 

Tbs  first  eight  coefficients  are  6^=  1.  666666  x  10*^, 

Oj  *  4.232804  x  10* 3.  e2  =  9. 472010  x  lO*4,  e}  *  1.  782967  x  10 ’4, 

•4  *  2.859895  x  10 ‘5,  eg  ■  3.978986  x  10‘6, 

•  4  *  4.878450  x  10"7.  and  e?  *  5.  341019  x  10“8. 

Tho  integrals  for  either  B^(T*)  or  JfJ(T*)  mast  be  done 

numerically.  Of  the  two.  the  numerical  integrations  for  B*  (T*)  are 

m  * 

mm%  easier;  the  numerical  details  are  described  In  Appendix  A. 

When  tbs  exponential  in  Eq.  (57)  is  expanded.  1*  (T  )  is  found  in 
terms  of  a  power  series  in  T*”1  with  the  coefficients  given  by 


46 


<  .  tjf  ('■'  K**  M(K*)  ctK* 

*  71!  4 

®- <T,)  “  t1*  T  *“*■  L"t^ 


Thus 


(81) 


(82) 


n*o 

-1 


=  -  4.8167  x  10'2  , 


We  find  that  iQ  =  2.  5400  x  10 

f,  *  8.  207ixl0‘3,  f„  *  -  1.  1886  x  10*3,  f.  =  1.4766  xl0“4, 

2  3  4 

f*  *  -  1. 5985  x  10"5,  f,  =  1.  5302  x  10’6,  and  f_  =  -  1.3115x10"! 
5  o  l 

Because  of  the  numerical  integrations,  it  is  uncertain  whether  or  not 

the  i  are  correct  to  five  significant  figures, 
n  17 

Recently,  Nosanow  and  Mayer  have  shown  that  the  reduced 
second  virial  coefficient  for  the  Lennard- Jones  (6-12)  potential  satis¬ 
fies  the  differential  equation 

t  *(S  (S3) 

Thu  equation  can  be  eolved  by  letting  B*  ■  T  $(T*  S  , 

*  *  -1 

and  noticing  that  £  (T  )  satisfies  the  differential  equation  of  a 

15  18 

confluent  hyp^r geometric  function  .  Thus 


e*(r*> -  t  (vr*-1*  f(Vu,v»./r»-')]  («■*> 

Here,  R  .  and  R^  are  constants  which  can  be  determined  by  com¬ 
paring  the  series  form  of  Eq.  (84)  with  a  result  previously  obtained 

19 

by  Lennard^Jeoes  : 


L.  H.  Nosanow  and  J.  E.  Mayer,  J.  Chem.  Phys.  28,  874  (1958). 
An  error  of  one  sign  in  the  differential  equation  given  by  Nosanow 
and  Mayor  has  been  corrected  in  £<§«  (83). 

18  * 

A%  alternative  method  of  arriving  at  £q.  (84)  la  given  by  L.  F. 
Epstein  and  Q,  M.  Roe,  J,  Chem.  Phys.  ^19,  1326  (1951), 

-Reference  3,  M>  T.  C.  L*»,  p.  163. 
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K,  -  -  2***  r  (-•/*«)  (85) 

rx  *  -  r,/m  ro/*»)  (86) 

Eq.  (84)  can  be  written  in  term*  of  the  Whittaker  functions^ 

M,  ,l/4(T*-1)  »nd  1/4<T*'*).  Thu*, 

B*fT»)  -  e<lT“)"[R1Ht.i(T-')  »RtMt  (87) 

By  using  the  expansion* 


M.  ,(T*-) 

i.'4 


r  £  _ ]■ 

L  j-o  I  ^o^fn-J)!  i  ji 


•*-* 


(88) 


and 


M  CT*'"; 

U 


-*~v* 


£  [i-t 

*»•  u  i«o  J 


2Vj  (SA)jCn-J)l 


7T] 


I*-* 


(89) 


4^ 

Eq.  (87)  become*  equivalent  to  the  result  given  for  B*(T  )  hy 
Nosanow  and  Mayer  except  that  they  give  recursion  relations  for  the 

♦  -  n 

coefficients  of  T  while  Eq*.  (88)  and  (89)  show  the  explicit 
form  of  these  coefficient*. 

The  reduced  *econd  virial  coefficient  for  this  potential  ha*  been 
20  *  * 

tabulated  .  Thu*.  B*jj(T  )  i*  found  by  using  Eq.  (3)  in  it*  reduced 
form.  The  contribution*  to  the  second  virial  coefficient  for  the 
Lennard- Jones  potential  are  shown  in  Table  IV  and  Fig.  10. 

With  the  table*  given  in  this  section,  it  i*  easy  to  calculate  the 
number  of  bound  and  metastably  bound  double  molecules  in  a  gas  if  the 
force  constanta  are  known.  If  the  activity  coefficients  are  set  equal 


a 


20 


References,  M.T.G.L.,  pp.  1114  and  1119. 


Contributlone  to  the  Second  Viriel  Coefficient 
for  th>  Lennerd-Jones  (6-12)  Potential 
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e 


I 


T* 


I 
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to  unity  in  Eq.  ( 14)  and  if  n^,  the  number  of  moles  of  single  mole¬ 
cules  is  replaced  by  n,  the  number  of  moles  that  would  be  present 
if  molecular  association  did  not  occur,  then  on  using  Eq.  (15)  the 
mole  fraction  of  double  molecules  is  given  approximately  by 

\  »  -tjB*  + 

while  the  mole  fractions  of  bound  and  metastably  bound  double  mole¬ 
cules  are  given  by 

Xlb-  -b.  6*  n/V 
Xw*  -faX»/v 

As  an  example  of  the  type  of  results  to  expect,  we  list  in  Table  V 
the  mole  fractions  for  a  number  of  gases  at  several  temperatures. 
The  calculations  are  for  the  Lennard-jdnes  potential  and  at  a 
concentration  of  1  mole/22.4  liters  (i.e.  n  =*  1,  V  =  2.  24  x  10* 


Table  V.  The  Hole  Fraction  of  Dimers  at  a  Total  Concentration  of  1  mole/22.4  liters 
Calculated  Using  the  Lennard-Jones  (6-12)  Potential. ^ 


(Continued) 


(Concluded) 
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V.  THE  LIFETIMES  OF  METASTAB LY  BOUND  MOLECULES: 

It  was  mentioned  in  the  introduction  that  the  mean  lifetime  of 
metastably  bound  molecules  will  have  an  effect  on  the  transport 
properties  of  gases.  This  is  because  double  molecules  with  a  life¬ 
time  less  than  the  mean  time  between  collisions  will  behave  more 
like  single  molecules  than  double  molecules.  It  will  be  seen  below 
that  except  at  low  temperatures  (the  comparison  is  being  made  at  a 
fixed  density)  most  metastable  molecules  have  lifetimes  sufficiently 
long  so  that  effects  of  dissociation  will  not  have  to  be  considered  in 
calculations  of  the  transport  properties. 

The  theory  of  the  mean  lifetime  of  metastably  bound  double 

21  22  23 

molecules  is  very  similar  to  the  one  body  model  theory  ’  '  of 

OC  -radioactivity.  There  are  several  ways  of  calculating  rasan  life¬ 
times  but  we  shall  employ  the  simplest  method  using  the  W.K. B. 
approximation.  This  is  the  quasi-ciassical  approach  discussed  in 
reference  22.  From  a  classical  mechanical  point  of  view,  the  meta- 
stable  molecule  has  a  vibration  period  t.  This  is  the  time,  as 
calculated  by  classical  mechanics,  for  the  system  to  move  once  back 
and  forth  across  the  potential  well.  Thus, 


t 


(90) 


H.  A,  Bethe,  Rev.  Mod.  Fhys.  9*  69  (1937),  See  especially  part 

XI. 

E.  C.  Kemble,  “The  Fundamental  Principles  of  Quantum  Mechan¬ 
ics"  (McGraw-Hill,  1937},  Section  31. 

D,  Bohm,  "Quantum  Theory  *  (Prentice -Hall,  1951)*  Ch.  12. 
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where  r  and  r_  are  the  two  turning  points  for  a  fixed  value  of  K 
1  2 

and  Kb^.  According  to  quantum  mechanic  a,  there  will  be  a  probability 
of  tranamiasio*V  that  a  particle  hitting  the  potential  barrier  at  r 2  will 
penetrate  the  barrier.  Thus,  the  probability  of  diasociation  in  one 
second,  P,  ia 


P  -  d/t 

and  the  nnean  lifetime  T  ia 


(91) 


X  -  i/p  -  t/e 

The  tranemiaaion  coefficient  T  can  be  found  by  the  W.K.B. 

.  .  .  .21 

and  ia  given  by 


(92) 

method 


e  - 

where 


(93) 


s  ’  *  r ,94> 

Here,  r  ia  the  turning  point  inaide  the  potential  well  (see  Eq.  (90)) 
and  r^  is  the  turning  point  outside  the  well.  For  the  W.  K. B.  method 
to  be  applicable,  G  should  be  somewhat  larger  than  one. 

The  mean  lifetimes  of  metaatably  bound  molecules  will  be  calcu¬ 
lated  for  systems  interacting  according  to  the  Lennard-Jones  (6-12) 

potential.  When  Eqs.  (78)  and  (7)  and  the  reduced  variables  intro  - 

*-2 

duced  at  the  beginning  of  part  IV  are  used  (with  y  =  r  ),  Eq.  (94) 
becomes 


6  *  jjr  J*4  -  iy*  7  K"b -  ** 

jfc 

Here,  .A  is  a  parameter  which  also  appears  in  the  quantum  mecbin- 

24 

ical  equation  of  state  of  a  Lennard-Jones  gas, 


24 


Reference  3,  M.  T.  G.  L.  ,  p.  421. 


(96) 


Alto,  with  thi*  change  of  variable*,  Eq.  (90)  for  t  become* 


The  right  hand  aide*  of  Eqa.  (95)  and  (97)  lead  to  hyper- 
elliptic  integral*  and  hence  are  evaluated  numerically.  The  detail* 

are  described  in  Appendix  B.  The  integral*  were  evaluated  for 

&  4*  4*2 

aevaral  value*  of  K  for  three  different  fixed  value*  of  K  b  . 

*  *2 

The  potential  curve*  corresponding  to  these  three  value*  of  K  b 

$  4 

are  shown  in  Fig.  2.  In  Table  VI,  2  AG  and  2  £  At/h  are  tabu¬ 
lated.  Thu*  if  the  potential  parameter*  of  a  ga*  are  known,  the  mean 
lifetime*  of  metaatably  bound  molecule*  can  be  eetimated.  Table  VI 
al*o  show*  the  mean  lifetime  for  metaatably  bound  atgon  molecule* 

(  A  *  0.  166,  <E/k  -  119-  6°K),  For  K*  -  0,49*  G  i*  approxi¬ 
mately  0.  5  in  the  caae  of  argon  so  that  the  W.K.  B.  approximation 

* 

doe*  not  apply*  but  for  tbe  other  value*  of  K  .  G  varies  between 

2  and  35  *o  that  the  W.K, B.  approximation  i*  valid.  $|ace  tbe  mean 

-9  -10 

time  between  colUaioaa  in  of  tbe  order  of  10  or  10  seconds 
at  one  atmosphere  pressure  and  ordinary  tempo ratur**,  on*  *••*  that 
meat  metastable  argon  molecule*  will  not  dissociate  before  under - 
goiag  e  collision.  Tbe  only  exceptions  will  be  molecules  in  meta- 
s table  etetes  wbicb  happen  to  have  energies  lying  close  to  the  top  of 
the  hnmp  in  the  potential.  However*  at  low  temperatures,  when  the 
moan  time  between  colli* lone  incr*a*e*,  a  larger  somber  of  meta¬ 
stable  molecules  wilt  dissociate  before  a  collision  with  another 
molecule.  This  will  load  to  a  rather  email  correction  in  the  trans¬ 
port  coefficient*  of  the  gas. 


IK*-- 


Tabic  VI.  The  Function*  2_A*  G  and  2€  A  t/h  for  Calculating 

the  Moan  Lifetime  of  Metaetabla  Double  Molecules  Using 
the  Lennard-Jones  (6-12)  Potential 


*  *2 
Kb 

* 

K 

2A*G 

2C  A?t/h 

(eec.) 
argon'  7 

.02 

8.18 

3.75 

4.75  x  107 

.64135 

.04 

4.74 

4.28 

5.25  x  10“l 

.06 

2.64 

5.15 

8.30  x  10'6 

i 

.08 

1.17 

6.96 

3.98  x  10’9 

.06 

12.0 

1.83 

2.19  x  1016 

1.3261 

.12 

6.38 

2.04 

4.43  x  104 

.18 

3.92 

2.41 

3.78  x  10"3 

.24 

1.74 

3.14 

3.82  x  10-8 

.40 

2.23 

2.29 

3.94  x  10‘7 

>6435 

.43 

1.49 

2.56 

8.52  x  10'9 

.46 

.824 

3.05 

2.76  x  10’10 

.48 

.199 

4.53 
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1.42  x  10 

VI.  THE  NUMBER  OF  VIBRATIONAL  LEVELS  OF  A  DOUBLE 
MOLECULE: 

During  the  course  of  this  research,  an  equation  for  the  maxi¬ 
mum  number  of  vibrational  states  of  a  dimer  was  obtained  for  the  case 
of  iooIoc'uIou  late raetlt\g  with  a  LtoinartLJtmen  (6fl2)  potential,  This 
equation  will  be  derived  in  this  section.  The  results  obtained  here  will 
justify  the  use  of  classical  mechanics  instead  of  quantum  mechanics 
in  the  evaluation  of  the  partition  functions  in  section  III.  It  is  known 
that  the  methods  of  classical  statistical  mechanics  are  applicable 
provided  that  aE  «  kT,  where  AE  is  the  energy  difference 
between  two  neighboring  quantum  states  of  a  system.  In  calculating 
^E  for  a  dimer,  only  the  change  in  vibrational  energy  need  be 
considered  because  the  rotational  states  are  very  close  to  each  other 

except  for  hydrogen  and  helium  isotopes.  Thus,  in  the  case  of  the 

$ 

Lennard- Jones  potential,  one  must  show  that  AE/e  «  T  .  For 

example,  it  will  be  seen  that  for  argon,  there  are  nine  vibrational 

levels  in  a  potential  well  of  depth  6  ,  so  that  classical  mechanical 

★ 

methods  are  valid  if  1/9  «  T  .  Except  for  hydrogen,  helium,  and 
neon,  the  results  below  show  that  quantum  mechanical  corrections 

sfc  % 

should  be  small  for  most  gases  at  values  of  T  as  low  as  T  =0.  5. 

* 

For  values  of  T  a  little  larger  than  one,  the  quantum  corrections 
for  neon  also  become  small. 

25 

According  to  the  W.K.B.  method  ,  the  vibrational  quantum 
number,  v,  is  given  by  the  equation 


25 


L.  I.  Schiff,  "Quantum  Mechanics",  (McGraw-Hill,  1955), 
Ch.  VII. 
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(V  *i)  -  ^7  (e  -  -  *££*£}'*  *  (98) 

The  limits  of  integration  are  the  two  positive  zeros  of  the  integrand, 

E  is  the  energy  eigenvalue,  and  X  is  the  angular  momentum  quan- 

-2  26  „  ,2 

turn  number.  The  coefficient  of  r  contains  (  /  +  2)  instead 

of  Jt{JL  +  1)  because  of  the  singularity  in  the  potential  at  r  =  0. 

Eq.  (98)  can  be  rewritten  in  terms  of  the  quantum  mechanical  para- 
❖ 

meter  J\_  ,  which  was  introduced  in  Eq.  (96),  and  the  reduced 

*  *  ,* 

quantities  r  and  <p  .  For  the  Lennard-Jonee  potential,  p  is 

given  by  Eq.  (78)  and 


-A "U*it 

lb  ir*-  r*x 


(99) 


The  integral  on  the  right  side  of  Eq.  (99)  is  hyperelliptic  and 
thus  the  energy  eigenvalues  cannot  be  obtained  except  by  numerical 
methods.  However,  if  we  are  only  interested  in  the  maximum  number 
of  vibrational  states,  the  integrand  can  be  simplified.  First,  the 
maximum  number  of  vibrational  states  occurs  when  =  0.  This  is 
because  the  addition  of  the  centrifugal  term  makes  the  effective  poten¬ 
tial  well  narrower  and  less  deep,  and  both  of  these  factors  decrease 
the  number  of  energy  levels.  Therefore,  we  need  only  consider  the 
case  where  JL  =  0.  Under  these  circumstances,  the  coefficient  of 


*-2  .  A* 

r  is  -=±- 


-2 

which  is  never  more  than  1.  50  x  10  (the  value  for 


69  7 T* 

3  -4 

He  )  and  for  most  gases  is  of  the  order  of  10  .  Thus,  only  a  small 


’  *1 

P.  M.  Morse  and  H.  Feshbach,  "Methods  of  Theoretical  Physics  , 
(McGraw-Hill,  1953)  Part  II,  p.  1101. 
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error  is  made  by  neglecting  the  term  containing  r  .  Secondly,  for 
a  bound  state,  E  cannot  be  greater  than  zero.  If  1?  is  set  equal  to 
zero  in  Eq.  (99),  the  right  side  will  in  general  no  longer  be  equal  to 
v  +  $,  but  the  greatest  integer  contained  in  the  right  side  of  Eq.  (99) 
minus  |  will  be  the  maximum  value  of  the  vibrational  quantum  number. 


Thus 


-i] 


(100) 


In  Eq.  (100),  the  integration  limits  ate  the  positive  zeros  of  the  inte¬ 
grand.  In  this  section  only,  the  square  brackets,  [  j  ,  indicate  the 
largest  integer  Less  than  or  equal  to  the  number  inclosed  by  the 
brackets. 

The  integral  in  Eq.  (100)  is  easily  evaluated  in  terms  of  the 

incomplete  elliptic  integral  of  the  first  kind^by  making  the  substitu* 

*.2 

tlon  y  m  r  .  Thus, 


-t] 


(101) 


where 


r  -  w'ffVrt) 

a  -  a  *  $  r 


(102) 


(103) 


The  maaimism  number  of  vibrational  states,  M,  is  eqaal 
(because  v  «  0  is  the  first  vibrational  stats).  By  inserting  the 
numerical  value  of  F(  V ,  k)  in  Eq.  (101),  wo  find 


H  *  [•**<*•  *i] 


(104) 


P.  T.  Byrd  and  hft.  0.  Friedman,  "Handbook  of  ilHytk  integrals 
tor  Engineers  and  Physicists"  (Berlin:  J.  tpflngsr,  1934). 
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Table  VII  shows  the  results  for  several  molecules.  It  is 

4 

interesting  to  note  that  for  He  these  calculations  predict  one  vibra¬ 
tional  level  lying  near  the  top  of  the  potential  well.  More  accurate 

methods  show  that  probably  no  discrete  levels  exist.  However,  this 
4 

result  for  He  indicates  that  in  the  case  of  the  JLennard- Jones  (6-  12) 
potential,  the  W.K.B.  method  does  give  a  fairly  good  approximation 
for  the  number  of  levels  even  when  the  maximum  quantum  number  is 
low,  while  for  most  potentials  one  can  only  use  the  W.K.B.  method 
with  confidence  when  dealing  with  high  quantum  numbers.  It  should 
be  noted  that  for  the  polyatomic  molecules  in  Table  VII,  the  number 
of  vibrational  levels  refers  to  vibrations  of  the  bound  molecules  and 
not  to  the  vibrational  levels  of  the  individual  molecules.  In  such  cases, 
it  is  assumed  that  there  is  no  coupling  between  the  vibrations  of  the 
molecules  and  the  vibrations  within  a  molecule. 


Table  VII.  The  Maximum  Number  of  Vibrational  Levels  for  Several 
Dimers  Using  the  Lennard-Jones  (6-12)  Potential 


Single 

Molecule 

A* 

1.684  1 

a*  2 

Maximum  Number  of 
Vibrational  Levels 
in  Dimer 

_  * 

He 

2.67 

1.13 

1 

Me 

.593 

3.34 

3 

A 

.186 

9.54 

9 

Kr 

.102 

17.0 

16 

Xe 

.064 

26.8 

26 

ch4 

.239 

7.54 

H2 

.226 

7.94 

7 

CO 

.220 

8.15 

8 

°2 

.201 

8.87 

8 

The  values  of  the  parameter  A  are  fro*  Reference  3 
M.T.G.L.,  d.  423. 
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APPENDIX  A.  THE  COMPUTATION.  OF  THE  COEFFICIENTS  fR 

OF  EQ.  (81)  FOR  THE  LENNARD- JONES  (6-12)  POTENTIAL: 

In  this  appendix,  the  numerical  calculations  used  in  finding  the 

coefficients  f  of  Eq.  (81)  are  described.  For  the  Lennard- Jones 
n 

(6-12)  potential,  Eq.  (58)  becomes 


A(*">rK»  i  f  _s_  _  i'V* 
ii(K  )  *  J  J  ' 

rfW 


(A  - 1) 


One  finds  K  bf  as  a  function  of  rh  by  solving  d  <f) gff  /dr  =  0  for 

*  *2  *  *  *  * 

K  b^  .  Then,  r^  is  found  by  solving  P  (rj1)  =  K  .  The  results 


(A  -2) 


and 

-  4^[|  -  [a  ♦  -»’f]  (A- 3) 

>:<  * 

Numerical  methods  are  used  to  find  r  .  The  quantity,  K  ,  which  is 

I.  c 

the  upper  limit  for  the  integrals  f^,  is  needed.  It  is  found  that 

jjc  t  if  2**2 

K  =0.8  by  solving  d  ^l^^/dr  ~  ®  and  ^  f>eU  /dr  =  ° 

simultaneously. 

The  integral,  M(K  ),  can  be  written  in  a  more  convenient  form 

*_2 

for  numerical  integration  by  making  the  change  in  variables  y  =  r 
Thus, 


±i 


/J.  '  /  i/  ^  c  ,  a' 

moo®  yJ  ( y-  -  Hy  )  - 


(A-4) 


and  in  this  form  it  is  seen  that  M(K  )  is  a  hyperelliptic  integral  in 


the  case  of  this  potential.  The  function  which  is  actually  calculated  is 
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¥  * 

M(K  j  -  the  reason  for  this  being  that  as  K  — *0,  M(K  )  has  a 

-I*  K* 

* 

logarithmic  infinity.  To  show  this,  we  shall  write  M(K  )  in  still 

•r_  3 

another  form,  Let  t  -  r  in  Eq.  (A-l). 


M(K  )  -  y  J  0  t  +■  Tt1  -  utH/s  )  t 


(A-5) 


$  * 

An  upper  bound  for  M(K  )  as  K  — ►O  is  found  by  noticing  that 


>  a!  -  Aibr 

-  h  t  *•  m*'* 


(A-6) 


The  equal  sign  holds  at  t  -  t  -  otherwise  Eq.  (A- 6)  is  an  inequality. 


Hence 


m(k*)  —  \  f*  (i  -tx  +  ~ 


(A-7) 


From  Eq.  (A -2),  it  follows  that  as  K  approaches  zero,  t  approaches 

*  I  •». 

(K  / 8)  .  \'so,  as  K  approaches  zero,  t^.  approaches  unity  because 

‘Jf. 

lim  if  (K  )  -  1.  The  right  side  of  Eq.  (A-7)  is  evaluated,  in  the  limit 

k'  — *0  7 

-t  2 

of  small  K  ,  bv  letting  x  ~  t  ,  integrating  by  parts  (the  division  being 

— ~  and  the  radical),  using  Eqs,  <236.00)  and  (236.  16)  of  Byrd 

X  S/J- 

27  t 

and  Friedman  ,  and  then  taking  the  limit  as  K  — e0.  Thus, 


1 1  to  M(K*)  —  -  -r-  |ti  K*  -  A-  ■+■  y  1 7)  2.  * 

$  * 

A  lower  bound  for  MfK  )  as  K  — >0  is  found  by  noticing  that 


(A -8) 


for  small  K  , 
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%/% 

K*  bjT*  _  JsT  *  CK^bTM  (A-9) 

4t"'*  Ht1-  ~  6(5K^)l/J-t 

/  •x  V  *  ,  , 

The  equal  sign  holds  in  Eq.  (A-9)  for  t  =  (  ;  otherwise  the 

inequality  holds.  (This  can  be  verified  by  the  usual  methods  of 
calculus.)  Therefore 


M(K*)  >  %  j  [*ta  -+t 


fK%r^  “i,/a‘  dt 

4  (  3  K*)  J  ta/a* 


(A -10) 


In  Eq.  (A -10),  the  limits  of  the  integral  are  the  positive  zeros  of  the 

radical.  These  zeros  lie  between  t,^  and  tf  because  tfa  and  tf  are 

the  zeros  of  the  radical  in  Eq.  (A-5)  while  the  integration  limits  of 

Eq.  (A -10)  are  the  zeros  of  a  radical  which  is  always  less  than  or 

equal  to  the  radical  of  Eq.  (A-5)  in  the  region  of  interest.  When  the 

right  side  of  the  Eq.  (A-10)  is  integrated  by  parts  (the  division  being 

dt  and  the  radical)  and  Eqs.  (256.00),  (256.  11),  and  (256.  12) 

t*/1L  27 

of  Byrd  and  Friedman  are  used,  we  find  on  taking  the  limit  as 

* 

K  — »-0  that 


1 1 fn  M(K*)  £  -  InK*  -  2  +■  3  )n  2.  (A-ll) 

K*-*0  ^ 

* 

Unfortunately,  a  better  estimate  of  the  limit,  lim  M(K  ),  could 
not  be  obtained  than  those  of  Eqs.  (A-8)  and  (A-ll).  Therefore,  in 
finding  the  value  of 

f  «  (ill*  fMV1’ M(K*)  dK*  (A  - 1 2) 

n  hi  J0 

* 

numerically,  the  technique  of  subtracting  the  limiting  form  of  M(K  ) 
from  the  integrand  of  Eq.  (A -12)  and  then  integrating  the  logarithmic 
term  analytically  cannot  be  used.  The  method  used  consists  in  writing 
Eq.  (A -12)  as 
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(z!)n 

n« 


K>n~  MfK*) 

-h  K - 


(’hK*1 


d  K* 


th?  (Xt  ^ 

M  4  -  /  ?!  K 


<A-13) 


whe  re 

X  =  -  K'  U  K*"  ‘  (A*14) 

When  n  ~  Or  the  integrand  of  Eq.  «A-13»  has  an  infinite  slope  at 

x  -  0.  However,  since  ^ —  -i-  ,  and  is  thus  finite  when  x  =  0, 

-  In  K  w  o 

this  does  not  cause  much  error  since  the  integrand  is  evaluated  for 

closely  spaced  values  of  x.  In  fact,  it  is  estimated  that  the  error 

-4 

from  this  source  in  f^  :s  somewhat  less  than  10  .  For  other 

values  of  n ,  this  difficulty  does  not  occur. 


The  integrand  of  Eq.,'A-13)  is  evaluated  at  0,  f  3-Xu_-  s 

Z^Xr,  »7Xc  m  Xc 

°  *  5  30  '  lio  '  I.V.’ 


a 

■  ,  3  ^ 

8  Xc. 

24  X, 

’  *  IO 

’  1  o 

'  5  *  *  ’  •’  to  ’ 

30 

5  30 

\ 

.  For 

each  of  tr.ese 

values 

Of  X; 

and  then 

-In 

was  found  fro 

m  Eq, 

;'A-4n 

was  computed  at  V,  ♦  n'v.  -  \'t  for  n  equal  to  0,-! — , -2— ,  ... 

n  I  •  h  60  60 

>8 


Z 

Zo 


20 


,  _i5jS_  ,  $6  .  .  .  s  1,  Over  the  intervals  n  --e  0  to 

20  60  60 


n  -  and  n  •  —  to  n 
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1  a  seven  point  integration  formula  can  be 
used  and  over  the  intervals  n  -  -rf;  to  n  =  and  n  "12.  to  n  ~  19u 

*- u  £0  OO  o  o 


IS 

a  nine  point  integration  formula**  can  be  used.  After  was 

“In  K* 

computed,  the  integral  of  Eq.  (A  -13}  was  found  by  using  a  seven  point 
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•'Tables  of  Lagrangian  Interpolation  Coefficients”,  Mathematical 
Tables  Project,  Federal  Works  Agency,  (Columbia  University 
Press  1948}.  n,  XXXII. 
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integration  formula  over  the  intervals  £o}  and  [  }  j  , 

a  six  point  formula  over  the  interval  £~fc  }  j  ,  and  a  five  point 

f  2.9  X  1 

formula  over  the  interval  [  ~3o"t  >  Xt  j  .  Table  A-l  shows  some  of 

the  intermediate  results  in  the  calculation  of  f  . 

n 

The  accuracy  of  the  numerical  integrations  is  difficult  to  deter¬ 
mine.  The  five  figures  given  for  f  may  not  all  be  significant. 
However,  at  worst,  the  error  should  not  affect  more  than  the  fourth 

figure  in  the  first  few  f  and  the  third  figure  in  the  last  few  f  . 

n  n 


APPENDIX  B.  THE  COMPUTATION  OF  THE  INTEGRALS  li  ! 

EQS.  (95)  AND  (97): 

The  integral  o£  Eq.  (95)  is  easily  evaluated  by  numeric  l  int-r- 

*  *  *2 

gration.  For  a  given  value  of  K  and  K  b  ,  the  integration 

limits  and  are  first  computed.  These  integration  limits  are 

the  two  smallest  positive  roots  of  the  radical  in  the  integrand.  The 

n 

integrand  was  computed  at  +  -g-  ( y^  *  v^)  for  n  equal  to  0,  ,  5,  1, 

r-v.  7 

1.  5,  2,  3,  4,  5,  and  6  for  all  values  of  K  and  Kb’  in  Table,'  V.  For 
*  ■<#  ’  Z 

K  equal  to  .02  and  .04  (Kb  0.64355)  and  K  equal  t<~>  ,06  and 
*  *2 

0.  12  (K  b  --  1.  3261)  the  integrand  was  also  computed  for  v  =  0.25 


and  0.  75.  The  integral  was  evaluated  by  using  a  five-pcmt  inegra- 
28 

tion  formula  over  the  intervals  n  =  0  to  n  2  and  n  -  2  to  n  6 
except  when  the  integrand  was  evaluated  at  i  ~  0.  25  and  n  ss  0.  75. 
In  the  latter  case,  a  five -point  integration  fo  rmula  was  used  over 
the  interval  n  ~  0  to  n  •-  1,  a  th:  ee-  point  formula  from  n  =  1  to 
n  2,  and  a  five-point  formula  from  n  s*  2  to  r»  -  6. 

To  evaluate  the  integral  of  Eq.  (97),  the  integration  limits  .ire 
computed  as  above.  The  limit  y .  is  the  largest  positive  root  of  the 
radical  in  the  integrand.  The  integrand  is  infinite  at  the  two  end¬ 
points  so  near  the  endpoints  the  polynomial  m  the  radical  ia  appro  <i- 

mated  by  (y  -  VqWYq)  .  Here,  f  (y)  ~  (K  *  -  4y°  +■  4y3  -  K*b*2y)/ 

_  i. 

(y  -  yQ)  and  yQ  indicates  either  y1  or  y  .  Now,  [f(y)l“2,  near 

the  endpoint  y^i  or  [-f(y)J  2,  near  the  endpoim  y^,  can  be 

expanded  in  a  Taylor  series  about  the  endpoint.  '1  he  apprcximatior 
(y  *  y0)f(y0)  wa*  uscd  f°r  ,he  polynomial  in  the  radical  as  long  as 

the  second  term  in  the  Taylor  series  expansion  of  j  v'v)J  2  or 

.  _ 1  1 

[-f(y)J  2  was  less  than  1%  of  the  constant  term  Whi  n  this 


approximation  Ui  used,  ike  contribution  to  the  integral  from  the  region* 
near  the  endpoiats  are  easily  evaluated  analytically.  To  obtain  the 
total  %alce  of  the  integral,  the  remaining  portion  between  the  new  end- 
points,  y'  and  ,  is  evaluated  numerically.  The  integrand  was 

*w4*aWd  *i  <y£  -+  ^  *y  '  -  ?*l  ior  n  equal  to  r0.  .% \,  l.  %  \  4, 

4.5,  5*5,  end  6.  A  five-point  integration  formula  was  used  from 
a  =:  0  to  n  =  2,  a  three -point  formula  from  n  s  2  to  n  =>  4, 
five-point  formula  from  n  «  4  to  n  a  6. 
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Th'.j  distance  of  closest  approach 
tion  of  b  for  a  fixed  value  of  the 
energy  K  less  than  K  . 


r  as  a  func ■ 
m 

initial  kinetic 


